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灰色 系统 GM(1,1) 与 三 角 函 数 耦 合 模型 的 研究 
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摘 要 : 灰色 系统 模型 由 于 其 具有 简便 、 数 学 背景 清晰 等 特点 ， 在 数学 建 模 与 近似 计算 等 领 
域 都 有 着 十 分 广泛 的 应 用 ,灰色 系统 模型 的 原理 可 以 概括 为 利用 时 间 响 应 函数 来 模拟 数据 序 
列 的 演化 ， 而 其 中 GM(1,1) 模 型 作为 最 基础 、 最 能 体现 灰色 建 模 思 想 的 一 个 模型 ， 应 用 最 广 
泛 。 传 统 的 GM(1,1) 模 型 的 时 间 响 应 函数 是 以 自然 常数 e 的 指数 构造 的 函数 ， 也 正 是 由 于 这 
个 特点 ， 所 以 GM(1,1) 模 型 具有 较 强 的 局 限 性 ， 在 实际 应 用 中 人 们 往往 会 将 GM(1,1) 模 型 的 
时 间 响 应 函数 做 一 定 的 修正 。 其 中 , 最 重要 的 一 种 修正 方法 就 是 将 原本 的 时 间 响 应 函数 与 其 
他 能 够 描述 数据 序列 性 质 的 函数 进行 耦合 ， 例 如 对 于 具有 一 定 周期 性 或 准 周期 性 的 数据 序 
列 ， 就 可 将 原本 的 时 间 响 应 函数 与 三 角 画 数 耦 合 ,构成 GM(1,1)- 三 角 函 数 耦 合 模型 。 本 文 便 
深入 探讨 了 这 种 模型 的 可 行 性 , 提出 了 耦合 模型 时 间 响 应 函数 中 参数 的 估计 方法 ,以 及 误差 
分 析 。 
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Abstract: Grey system model is widely used in mathematical modeling and approximate calculation 
because of its simplicity and clear mathematical background. The principle of grey system model 
can be summarized as the use of time response function to simulate the evolution of data series. 
Among them, GM (1,1) model is the most basic and can reflect the idea of grey modeling most 
widely. The time response function of the traditional GM (1,1) model is constructed by the exponent 
of natural constant e. Also, for this characteristic, the GM (1,1) model has strong limitations. In 
practical application, people often modify the time response function of GM (1,1) model. One of 
the most important correction methods is to couple the original time response function with other 
functions that can describe the properties of data series. For example, for data series with certain 
periodicity or quasi periodicity, the original time response function can be coupled with 
trigonometric function to form GM (1,1) - trigonometric function coupling model. In this paper, the 
feasibility of this model is deeply discussed, and the estimation method of the parameters in the time 


response function of the coupled model and the error analysis are proposed. 
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了 路 


1 引 


自 80 年 代 邓 聚 龙 教授 首次 提出 灰色 系统 理论 以 来 ， 灰 色 系 统 理论 在 各 大 数学 建 模 和 控 
制 论 领 域 都 逐渐 占有 一 席 之 地 。 它 不 同 于 大 数据 这 种 处 理 大 量 数据 的 学 科 , 而 其 诞生 就 伴 有 
“对 少量 数据 的 模拟 ”的 目标 背景 上 。 所 以 在 掌握 信息 少 ， 信 息 相 对 孤立 的 情况 下 ， 灰 色 系 
统 理论 成 为 了 建 模 时 出 类 拔 萃 的 选择 。 
但 是 ， 由 于 应 用 最 广泛 的 GM 模型 的 数学 结构 大 量 采 用 了 指数 函数 ， 导 致 该 类 灰色 模 
型 对 振荡 序列 的 模拟 始终 不 尽 人 意外 。 这 主要 是 源 于 指数 函数 的 函数 性 质 带 来 很 多 限制 。 对 
此 人 们 采取 了 很 多 办 法 , 例如 引入 缓冲 算 子 等 数学 技术 钻 ,， 但 这 些 都 不 能 从 根本 上 解决 函数 
单调 性 带 来 的 建 模 限制 。 
后 来 , 人 们 发 展 出 将 传统 灰色 模型 所 用 到 的 函数 与 工程 数学 中 的 一 些 函 数 耦合 起 来 的 技 
术 ， 这 起 到 了 意 想不到 的 效果 人 +。 这 是 一 条 弥补 传统 灰色 系统 模型 缺陷 的 可 行 的 路 ， 本 文 
也 是 基于 这 种 思路 将 传统 灰色 系统 模型 和 周期 性 函数 结合 起 来 , 用 于 模拟 具有 周期 性 和 准 周 
期 性 的 序列 。 
T= 为 简便 起 见 ， 本 文 就 以 最 常用 又 最 基础 的 GM(1,1) 模 型 与 三 角 函 数 进 行 耦 合 ， 所 得 到 的 
> 模型 称 为 GM(1,1)- 三 角 函 数 耦 合 模型 。 本 文 首 先 根据 灰色 系统 理论 给 出 其 时 间 啊 应 函数 , 再 
利用 一 些 数学 知识 给 出 了 模型 的 参数 估计 技术 。 并 结合 一 个 应 用 实例 对 有 具体 细节 给 予 了 说 
明 。 最 后 对 模型 进行 了 讨论 。 


2 GM(1,1)- 三 角 函 数 耦 合 模型 


传统 的 GM(1,1) 模 型 是 将 原始 序列 Xo ={ziza, zs 进行 1-AGO 生成 操作 ， 构 造 其 1- 
AGO 序列 了 。 之 后 利用 一 阶 微分 方程 作为 白化 方程 求解 时 间 响 应 函数 ， 得 到 对 Xi 序列 的 
"~ 模拟 ， 再 还 原 到 对 Xo 的 模拟 ， 从 而 实现 建 模 与 误差 分 析 。 但 正 是 因为 其 采用 线性 的 一 阶 微 
PP 分 方程 作为 其 白化 方程 ， 才 导致 GM(1,1) 模 型 的 时 间 响 应 函数 是 e4@9 型 的 ， 这 种 全 局 单调 性 
的 函数 严重 影响 了 GM(1,1) 模 型 的 适用 范围 , 对 于 震荡 序列 和 摆动 序列 的 模拟 显得 无 能 为 力 。 
对 于 震荡 序列 ， 已 经 发 展 出 了 基于 灰 数 理论 的 比较 普遍 的 方法 四， 所 以 本 文 重点 讨论 对 于 周 
期 或 准 周 期 摆动 序列 的 修正 模型 。 


2.1 时 间 响 应 函数 


我 们 直接 对 GM(1,1) 的 白化 方程 做 出 如 下 修正 
Xt) ta Xt) =Xsin(wt tp) (1) 


也 就 是 将 原来 等 式 右 边 的 灰 作 用 量 5 蔡 换 成 了 三 角 函 数 。 从 微分 方程 的 角度 来 讲 ， 这 样 做 的 
意义 在 于 把 原来 解 空间 平 直 的 超 平 面 构造 成 弯曲 的 , 使 得 方程 的 解 函 数 具有 周期 性 。 由 此 也 
引出 我 们 对 于 周期 性 和 准 周 期 性 的 定义 : 周期 性 序列 是 指 (1) 中 入 取 定 值 的 序列 ,而 准 周期 序 
列 是 指 (1) 中 入 时 间 依 赖 的 序列 。 下 面 以 三 幅 示 意图 说 明 一 下 常见 的 三 种 周期 性 或 准 周期 性 序 
列 。 
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a. 周期 性 序列 1 b. 周期 性 序列 2 c. 准 周期 性 序列 


图 1， 周 期 性 与 准 周期 性 序列 示意 图 从 图 中 来 看 ， 周 期 性 与 准 周期 性 序列 比较 直观 
的 区 别 在 于 ， 周 期 性 序列 的 周期 是 稳定 的 ， 而 准 周期 序列 的 周期 是 随时 间 变 化 的 。 
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那么 周期 性 序列 在 注意 初 值 条 件 X(t)| -= 和 (1) 的 情况 下 ， 微 分 方程 0) 有 通 解 
Xi(t)=Ae “+A(asin(wt+i+w)— weos(wtt wy)) 
A=[X(1)—AM(asin(w+ wy)— wceos(w+t wo))le" 
对 于 准 周期 性 序列 ， 则 要 结合 二 入 (t) 有 具体 的 形式 求解 通 解 。 如 果 和 (t) 的 线性 的 ， 那 么 求解 
并 不 是 困难 的 事 , 与 结果 (2) 没 有 本 质 上 的 区 别 , 而 如 果 和 (t) 是非 线性 的 , 那么 就 需要 采取 微 
扰 近 似 的 思想 ， 利 用 格林 函数 进行 求解 ， 这 里 不 再 袭 述 [9。 
而 利用 灰 导 数 的 定义 ， 我 们 有 Xo(z) 三 Xi(z) 一 Xi(z 一 1)， 也 就 是 说 
Xo(zZ) = = (eeos [e+(s— $)o|+ osin|y 十 (2- 3)|) (3) 


O) 


这 里 一 共有 四 个 参数 {a,w,Pp, 和 } 。 从 数学 形式 上 来 看 ， 这 是 两 个 函数 耦合 的 问题 ， 即 (3) 的 两 
项 ， 第 一 项 为 指数 项 ， 第 二 项 为 周期 项 。 而 各 自 的 系数 参数 4 和 入 则 代表 着 两 个 函数 在 整体 
中 所 占 的 比重 , 当 4 = A(a,w, wp, 入 较 大 时 , 函数 的 单调 性 更 加 显著 , 而 入 较 大 时 函数 的 周期 
性 更 加 显著 。 而 剩 下 的 参数 w,yp 继承 了 三 角 函 数 中 的 数学 意义 ， 是 序列 的 震动 频率 与 相位 。 
为 了 方便 后 续 讨论 ， 我 们 将 (3) 形 式 化 地 记 为 

Xo(7)=f(7)+2Ng(7) (4) 
其 中 f(x) 代表 指 数 项 ， 决 定 着 整个 函数 在 z 域 上 的 发 展 趋势 。g(zx) 是 周期 项 ,决定 着 整个 函 
数 的 振荡 趋势 .而 2 入 是 两 个 项 的 耦合 系数 , 在 两 个 项 构成 函数 的 过 程 中 有 权重 的 数学 意义 。 


2.2 参数 估计 技术 


在 找到 了 GM(1,1)- 三 角 函 数 克 合 模 型 的 时 间 响 应 函数 后 ， 下 一 步 的 工作 便 是 对 其 参数 
的 确定 提供 方法 以 便 实际 应 用 。 
根据 (4) 中 两 个 部 分 的 函数 性 质 ， 周 期 项 g(x) 决 定 着 整个 时 间 响 应 函数 的 振荡 趋势 ， 可 
以 考虑 将 周期 或 准 周期 理解 为 三 个 连续 零点 的 区 间 长 度 。 那 么 在 一 个 比较 特殊 的 条 件 
a,p 二 0 下, g(x) 的 零点 有 如 下 表达 式 
1 


和 (5) 


这 样 ， 三 个 连续 零点 的 区 间 长 度 仍然 是 我 们 所 熟悉 的 周期 == 2x/w。 对 于 条 件 a,p 二 0， 可 
以 证 明 ，a 和 gy 对 整个 9(z) 的 影响 仅仅 在 于 图 像 的 左右 平移 ， 并 不 会 影响 函数 局 域 的 性 质 。 
所 以 ap 三 0 是 在 探究 周期 表达 式 时 可 以 成 立 的 条 件 。 

既然 如 此 ,对 于 给 定 的 一 组 数据 序列 , 可 以 根据 其 变化 规律 , 找到 周期 或 者 平均 周期 卫 ， 
再 利用 了 = 二 2/w 直接 对 ww 做 出 估计 。 

在 平均 周期 和 w 确 定 以 后 ， 考虑 到 yp 具有 相位 的 意义 ， 通过 判断 图 像 的 平移 情况 就 能 对 
参数 gp 做 出 估计 。 正 如 我 们 在 三 角 函 数 中 所 做 的 一 样 。 具体 的 操作 将 在 后 续 的 应 用 实例 中 给 
出 。 


之 后 便 是 参数 a 的 确定 。 由 于 a 在 灰色 系统 模型 中 是 发 展 系 数 ， 决 定 着 整个 函数 的 发 展 
趋势 ， 所 以 它 的 确定 可 以 按照 极 大 值 拟 合 的 方法 进行 。 也 就 是 ， 取 序列 的 各 个 波峰 值 进行 传 
统 的 灰色 系统 建 模 ,确定 发 展 系 数 。 至 于 传统 的 灰色 系统 模型 ,这 里 不 再 玖 述 ， 可 查阅 相关 
文献 11。 


最 后 一 个 参数 ， 也 是 克 合 模型 的 灵 瑰 所 在 ， 灯 合 参数 入 。 由 于 指数 模型 的 单调 性 是 无 法 
避免 的 问题 ， 在 一 些 反 常数 据 区 间 ( 即 不 遵循 指数 律 的 数据 区 间 ) ， 需 要 耦合 参数 入 让 三 角 函 
数 的 占 比 增 大 ， 从 而 弥补 指数 模型 单调 性 的 缺陷 。 而 它 的 确定 也 并 不 十 分 容易 。 在 我 们 的 工 
作 中 准备 采取 变 分 法 来 对 其 做 出 估计 。 

当 @a,w,p 都 确定 时 ，Xo(z) 成 为 仅 依赖 于 zx, 入 的 函数 ， 不 妨 记 为 F(z, 入 )。 那 么 Xo 作为 
灰色 系统 给 出 的 模拟 序列 ， 我 们 希望 它 与 原始 序列 X(z) 的 相对 误差 尽量 地 小 ， 这 可 以 归结 
为 一 个 数学 问题 。 假 设 原 始 数据 序列 有 个 数据 ， 定 义 相 对 误差 累积 函数 
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2 Xn) ® 
虽然 求 和 记号 更 加 符合 实际 情况 , 但 求 和 所 定义 的 函数 并 不 利于 我 们 处 理 某 些 问 题 。 所 以 也 
可 以 对 原始 数据 做 一 次 精度 非常 高 的 拟 合 后 ， 以 积分 的 形式 将 它 推 广 成 一 个 泛 函 
NF = 
ov=/ Ys 本 dz (7) 


此 时 再 应 用 变 分 原理 6Qoy/6z 二 0 可 以 给 出 合适 的 入 ， 具 体 的 操作 将 在 应 用 实例 中 给 出 


llo 


3 应 用 实例 


任 取 一 组 具有 准 周 期 性 的 数据 序列 如 下 


表 1 原始 数据 
2 4 7 6 3 5 8 11 9 6 10 14 12 8 13 16 19 17 14 


根据 局 域 极 大 值 之 间 数 据 间 隔 可 以 判断 ， 这 个 数据 序列 具有 平均 周期 T= 二 14/3 。 由 此 
定 出 w= 二 27/T= 37/7 。 为 方便 起 见 ， 我 们 先 对 这 组 序列 的 发 展 系 数 a 进行 确定 。 
我 们 算 选 出 几 个 极 大 值 进 行 灰 色 系 统 建 模 ， 即 对 
X(3)=7,X(8)=11,X(12)=14,X(7)=19 
进行 GM(1,1) 建 模 。 最 终结 果 如 下 表 所 示 


表 2 极 大 值 的 GM(1,1) 模 型 


原始 数据 模拟 值 相对 误差 
\ \ 
11 10.724 2.51% 
14 14.151 1.08% 
19 18.672 1.72% 


发 展 系 数 a 二 -0.277 


之 后 ,我 们 需要 对 wp 进行 确定 。 根据 三 角 函 数 的 知识 ， 数 据 序列 的 第 一 个 极 大 值 点 出 现 
的 位 置 , 在 g = 二 0 时 理论 上 应 当 是 x 二 2.33 (x 应 当 取 正 整 数 , 但 这 里 为 了 估计 参数 ， 所 以 先 
暂时 取 数 值 计 算 结 果 )， 而 实际 上 是 z= 二 3， 相 当 于 整个 函数 右 移 了 。 那 么 可 以 求解 出 
p=-0.67. 

最 后 是 耦合 系数 入 的 确定 。 我 们 绘制 (6) 式 所 定义 的 函数 


> 


入 
图 2. 相对 误差 累积 函数 与 耦合 系数 关系 示意 


pa 


我 们 求解 到 最 小 值 出 现 的 点 对 应 着 入 = - 0.837 。 自 此 ， 模 型 中 所 有 的 参数 都 已 确定 ， 我 们 
带 入 具体 的 序数 与 原始 数据 序列 做 对 比 


表 3 耦合 模型 与 原始 数据 比较 


原始 数据 2 4 7 6 3 县 so i 16 19 17 14 
= 模拟 值 \ 373 701 632 281 ... e1721 18.14 15.93 12.74 
> 相对 误差 \ 68% 01% 53% 63% ... Bh 76% 45% 63% 9.0% 
yp 平均 相对 误差 5.6% 


表 3 说 明 耦 合 模型 得 到 了 一 个 非常 好 的 结果 , 平均 相对 误差 在 一 个 比较 低 的 水 平 。 而 如 果 直 
接 采 用 传统 的 灰色 系统 模型 对 原始 序列 建 模 ， 平均 误差 将 达到 48.3%。 这 并 不 让 人 惊讶 ， 因 
为 传统 的 灰色 系统 模型 本 身 在 振荡 序列 的 预测 上 就 有 很 大 的 缺陷 ， 属 于 建 模 机 制 问题 。 

本 节 提 供 了 一 个 具体 的 例子 说 明 耦 合 模型 中 一 些 具体 的 细节 步 又， 特别 是 参数 的 确定 。 
这 样 发 展 出 来 的 模型 弥补 了 传统 的 灰色 系统 模型 对 于 振荡 序列 模拟 和 预测 的 不 足 , 对 于 本 文 
所 定义 的 周期 与 准 周期 序列 具有 十 分 有 效 的 建 模 机 制 。 


4 讨论 与 改进 


即便 本 文 所 建立 的 耦合 模型 得 到 了 一 定 程度 的 应 用 ， 但 其 所 存在 的 不 足 仍 然 是 很 明显 
PE 的 。 


Sg 


首先 是 在 耦合 系数 的 确定 上 , 采用 变 分 法 得 到 的 系数 是 一 个 定 值 , 这 会 让 本 文 建立 的 类 
合 模型 仅仅 起 到 了 一 种 “缓解 ”的 作用 ， 而 不 能 从 根本 上 避免 时 间 响 应 函数 的 单调 性 所 带 来 
的 模拟 限制 。 比 如 我 们 研究 发 现 ， 当 耦合 系数 取 定 值 时 ， 整 个 模拟 函数 的 发 展 趋势 仍然 是 单 
调 的 , 即便 在 局 域 上 具有 很 丰富 的 振荡 结构 。 但 是 当 耦 合 系数 取 成 一 个 随时 间 发 展 而 变化 的 
函数 时 ， 事 情 发 生 了 根本 性 的 变化 。 由 于 耦合 函数 具有 权重 的 意义 ,在 那些 违反 数据 序列 发 


展 总 态势 的 区 间 里 (例如 股票 整体 是 上 涨 的 , 但 存在 局 域 性 、 剧 烈 的 下 跌 )， 就 需要 一 个 很 大 
的 耦合 系数 并 配合 三 角 函 数 的 升降 性 来 模拟 出 这 种 区 间 , 同时 指数 函数 的 参与 可 以 保证 发 展 
总 态势 的 模拟 。 


其 次 , 本 文 建立 的 耦合 模型 还 可 以 发 展 出 更 加 丰富 的 推广 , 例如 可 以 引入 对 数 函 数 来 控 
制 变化 幅度 ,引入 多 项 式 函 数 来 使 得 模型 对 局 域 极 值 的 模拟 更 为 准确 等 等 。 这 些 都 是 能 够 实 
现 且 能 够 得 到 广泛 应 用 的 推广 模型 。 

数学 建 模 是 一 门 深奥 而 又 直观 的 学 科 , 它 与 我 们 的 日 常生 活 奶 息 相 关 , 依赖 着 复杂 而 有 
趣 的 数学 理论 ， 它 是 我 们 认识 自然 规律 和 社会 规律 的 一 种 方式 。 
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